Finite Range Effects in Energies and Recombination 
Rates of Three Identical Bosons 



P K S0rensen, D V Fedorov, A S Jensen and N T Zinner 

Department of Physics and Astronomy - Aarhus University, Ny Munkegade, 
bygn. 1520, DK-8000 Arhus C, Denmark 

Abstract. We investigate finite-range effects in systems with three identical 
bosons. We calculate recombination rates and bound state spectra using two 
different finite-range models that have been used recently to describe the physics 
of cold atomic gases near Feshbach resonances where the scattering length is 
large. The models are built on contact potentials which take into account finite 
range effects; one is a two-channel model and the other is an effective range 
expansion model implemented through the boundary condition on the three-body 
wave function when two of the particles are at the same point in space. We 
compare the results with the results of the ubiquitous single-parameter zero-range 
model where only the scattering length is taken into account. Both finite range 
models predict variations of the well-known geometric scaling factor 22.7 that 
arises in Efimov physics. The threshold value at negative scattering length for 
creation of a bound trimer moves to higher or lower values depending on the sign 
of the effective range compared to the location of the threshold for the single- 
parameter zero-range model. Large effective ranges, corresponding to narrow 
resonances, are needed for the reduction to have any noticeable effect. 
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1. Introduction 

Low-energy quantum mechanical bound states present a number of surprising results 
when the state contains three or more particles. The most famous case is the Efimov 
effect that occurs when a three-body system has two-body subsystems that have a two- 
body bound state with zero energy. In this situation one can mathematically prove 
that an infinite number of three-body states appear when the particles are bosons [1] . 
More generally, the effect can also take place for two-component fermionic systems 
whenever the mass ratios are large enough [31 [31 H]. In recent years, it has become 
clear that few-body effects such as this can be studied in great detail in experiments 
with ultracold atomic gases via tunable Feshbach resonances [3 [B] . Producing a two- 
body bound state at zero energy is thus routinely done and signatures of three- and 
even four-body resonances have been observed [7ll8ll9l[T0 l [TT | [T2 l [T3 l fT4 l [T5 l fT6 l fT7 ] 

The few-body features that are studied in ultracold atom experiments are often 
identified through the rate at which atoms are lost from the experimental trapping 
potential. In fact, the densities and lifetimes of typical Bose-Einstein-Condensates 
(BEC) are limited by loss effects, primarily due to three-body recombination, a process 
where three particles interact and create a bound system of two particles (dimer) and 
the third particle carries away excess energy and momentum, generally resulting in a 
loss of all three particles from the confining trap. The loss rate is given by n = —an^ 
where n is the particle density and the recombination coefficient is a = C{a)ha'^ /m 
with C{a) a log-periodic function of the two-body scattering length a with period 22.7 
[2] . A strongly interacting BEC can be created by use of Feshbach resonances where 
the scattering length can be tuned using an external magnetic field 6 . 

The a'^ scaling of the recombination coefficient is valid when the scattering length 
is much larger than the range of the inter-atomic potential. Here we use different 
models to estimate the effect of the finite range of the interactions. All the models 
we use are based on contact interactions but, in contrast to the typical one-parameter 
implementation, we also incorporate a non-zero effective range. This is done by either 
modifying the Bethe-Peierls boundary condition using the effective range expansion 
or by using a two-channel model which not only has an inherent effective range but 
also qualitatively describes the physics of Feshbach resonances. 

The hyperspherical adiabatic approach is used to handle the three-particle 
problem. This method works equally well for the three scattering models we consider 
(zero-range contact interaction, effective range expansion and two-channel contact 
interaction). In order to calculate the recombination rate, we use the WKB-method 
with hidden crossings [18^ that takes the three-body system from the initial in-coming 
scattering channel into the atom plus dimer channel of lower energy via an excursion 
into the complex plane. Since this approach depends critically on the presence of 
a bound dimer, it is only applicable to the positive a side of a Feshbach resonance 
and we only consider recombination rates for a > 0. Along the way, we discuss the 
hyperspherical potentials for the different models, and also the scaling properties of 
consecutive Efimov resonances when effective range corrections are included. 

In order to also address the negative a side of the resonance, we calculate 
the spectrum of bound low-energy trimers within different models and study their 
dependence on the effective range corrections. Recent experiments have shown that the 
lowest three-body Efimov state has a breakup threshold at the three-atom continuum 
on the a < side of the resonance that seems to be uniquely determined by the 
background van der Waals parameter of the atomic system under study [15] . In 
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fact, the threshold, a^—), is experimentaUy found to fulfil a'^—) ^ — 9.8rvdWj where 
rvdw is the two-body van der Waals length scale [T^. This finding has induced a 
number of theoretical studies that aim to explain the proportionality and calculate 
the constant factor from basic knowledge about the two-body atomic potentials 
[m [20l [22l [211 US [21] ■ It can be most easily understood as consequence of the 
hard-core repulsion of the atom-atom interaction (typically of Lennard- Jones type 
with a corresponding van der Waals tail). This fact induces a hard-core repulsion 
in the three-body potential around r^dw as shown numerically in Ref. }21j and soon 
after by analytical means in Ref. [35] • Here we calculate the three-body spectrum 
for different scattering models and study the influence of effective range corrections 
on the threshold and the scaling properties of the states. Finite range corrections 
beyond the scattering length approximation have been considered before within various 
approaches [25l|26l|2a[Ml[2ai30l|3ll[32l[33ll34l|35]. However, to the best of our 
knowledge systematic studies using simple zero-range models have not been presented. 
Since zero-range models are extremely convenient in both few- and many-body studies, 
it is important to gauge their applicability which is one goal of the current study. 

The paper is organized as follows: In Sec. [2] we first present the hyperspherical 
adiabatic method which constitutes our theoretical framework and discuss the 
differences in the solutions to the hyperspherical three-body potential within the 
different two-body interaction models that we employ. Section [3] discusses the 
recombination rate within the different models using the hidden crossing technique. 
In section [4] we discuss the spectrum of three-body bound states with the various 
models and map out the dependence on the effective range. Section [5] contains a short 
summary along with conclusions and outlook. 



2. Formalism 



We consider a system of three identical bosons of mass m using hyperspherical 
coordinates defined from the Cartesian coordinates ri,rj,rk, of particles i,j, k as 

/2 /^^^ _ ^ 



^/2 ' V 3 V 2 

p2 = |a;,|2 + |y,|2, tana, = {^, (2) 

where fc} are cyclic permutations of {1,2,3}, p is the hyperradius and ai is a 
hyperangle [2|. The directions of Xi and yi comprise four additional hyperangles 
which, along with ai, are denoted collectively as f2. The hyperradius p is independent 
of the choice of {i, j, fc}. 

The wave function is expanded on adiabatic basis states <l'„(p, fJ). In the 
hyperspherical adiabatic approximation only the first term in this expansion is kept, 
this has been proven to be a good approximation when dealing with Efimov three-body 
states (trimers) [2]. The wave function is then 

vI'(p,r!) = p-5/Vo(p)$o(p,f^), (3) 
where $o(/0; ^) is the solution to the hyperangular equation 

^ + '^o^^' ^) = Mp)Mp. (4) 



Finite Range Effects for Three Bosons 



4 



where A is the grand angular momentum operator in hyperradial coordinates (see [5]) 
and Aq is the corresponding eigenvalue. The hyperradial function /o(p) is a solution 
to the hyperradial equation 

/ Ao(p) + 15/4 2mE\ „ , . „ ... 

7 Qoo(p) - -^j /o(p) = , (5) 

with 

$0) , (6) 
/ n 

where brackets indicate integration over all hyperangles. Numerically Qqq is found 
to be extremely small compared to other terms and we will not include it in further 
calculations. In the following we will suppress the subscript 0. 



Qo 



$0 



2.1. Two-body Potential Models 

Our main concern is the introduction of finite range effects beyond the single- 
parameter zero-range approximation (where only the scattering length is included) 
hut still using only contact interactions. There are different ways to do this and here 
we compare the results obtained by using a boundary condition or range expansion 
model and a two-channel model, and compare to the usual zero-range approximation 
with only the scattering length. To make our discussion self-contained and well-defined 
we now proceed to introduce first the scattering length only zero-rang model and then 
the two multi-parameter models. 



2.1.1. Zero-range Model Zero-range models generally use a contact interaction 
potential which is defined by a boundary condition on the logarithmic derivative of 
the wave function at zero separation. In hyperspherical coordinates this boundary 
condition becomes [23 



a(a,$) 



dai 



Qi=0 



(7) 



where a,; is the scattering length between particles j and k. Since all particles are 
equal we will suppress indices on scattering variables such as a. Using Faddeev 
decomposition of the angular wave function with s-states only, 



where 

is a solution to dH), 
eigenvalue equation 



N{p) sin [v - I ) 



(8) 



(9) 



= A + 4 and A^(/3) is a normalization factor, we find the 



%/3' 



(10) 



after rotating two of the Faddeev components into the coordinate system of the third 
pS] . For large positive pja there is an asymptotic solution of the form v i\plpja 
yielding a dimer binding energy of 
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(in units oi h — m = 1). For negative a there are no bound dimers. The hmit 
p/a = yields the solution v = is^ with sq = 1.00624. These are the basic aspects 
of the simplest single-parameter zero-range model when applied to a system of three 
identical bosons. 



2.1.2. Effective Range Expansion The first method for including finite range effects 
beyond the scattering length is the effective range expansion (range expansion or 
boundary condition model for short). It is a generalization of the boundary condition 
([7]). From scattering theory the effective range expansion of the phase shift 5 as a 
function of wave-number k is 



lim fccot(5(fc) = - - + -Rk^ 
k^Q a 2 



(12) 



where a is the scattering length and R is known as the effective range. This modifies 
the boundary condition ([7]) to (see [5S]) 



9(a^$) 



dai 



V2p 



Qi=0 



2^ V 



(13) 



where Ri is the effective range between particles j and fc, and the momentum in ()12|) is 
given by fc = Jy/iVip) [25]. A gain assuming that all particles are equal, the eigenvalue 
equation ([T0|) becomes 



V2p 



1 ^j^f" 
a 2 \y/2p 



(14) 



Inclusion of the effective range yields the dimer binding energy 

2 




(15) 



for \R\ <C a. Thus the dimer system can be more or less bound depending on the sign 
of the effective range. In the case of atomic Feshbach resonances the effective range is 
negative [6l, yielding less bound dimers. 

This model breaks down for sufficiently large positive effective ranges since there 
exist no solution to the eigenvalue equation (fT4|) for p < i? if i? is positive. To remedy 
this deficiency an additional known as the shape parameter, P, can be included in 

unD [55] 

liiR k cot S{k) = -- + l-Rk^ + PR^k^ , (16) 



and (|T3)) gets a similar additional term. Typical values of P are around 0.1. However, 
we note that results from binding energy calculations are largely insensitive to the 
exact value of P [25]. We therefore fix P = 0.1 in this study. 



2.1.3. Two- channel Model The other model for including finite range effects that we 
consider is a two-channel contact interaction model that takes the internal degrees 
of freedom of the atoms into account. This model has background scattering lengths 
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in open and closed channels that we denote aopon a-nd Odoscd respectively. The full 
scattering length and effective range parameters of this model are [36] 

^ "'opcn ~ 

^ ^closed 

R = — — : (18) 



k(k 



where /3 parametrises the coupling between the channels and k is given by the energy 
separation E* — h'^K^/2m* between the channels, where m* = ^ is the two-body 
reduced mass. The expressions are valid when a ^ (equivalently E* ^ fi^ /2m*a^ 
), which is always the case near a Feshbach resonance. At the resonance, i.e. a — oo, 
the effective range is [6] 

where AB is the width of the resonance, A/i the difference in magnetic moments 
of the channels and abg the background scattering length (the value away from the 
given Feshbach resonance). The resonance width and position Bq as given in the 
phenomenological expression [B] 

a{B) « (^1 - -^^^ , (20) 

which can be derived in the two-channel model. The parameters aoponj '^cioscdj /3 and 
E* can now be replaced by AB, Bq, the difference in magnetic momenta A/i and off- 
resonance scattering length ai,g to relate all quantities to physically measured values. 
Since we address neutral atom systems here the long-range atom-atom interaction is 
of the van der Waals type and can be characterized by the van der Waals length, rvdw- 
It has been shown that the scattering length in a potential with a van der Waals tail 
is ~ 0.956rvdw [37]. For simplicity we will assume that atg = rvdw and use the two 
interchangeably. This also means that we are working with abg > throughout this 
paper. We have checked that the case with at,g < gives the same qualitative results. 
Off resonance the model predicts the effective range R in terms of Rq and a 

Ria)=Ro[l-^y , (21) 

similar to the results obtained in [351 IMl HOI SD HlI Note that the effective range is 
always negative since Rq is always negative for an atomic Feshbach resonance W ■ Since 
the a-dependence of R of ((2T|) should be generic for Feshbach resonances (independent 
of the particular scattering model used), we will use R = R{a) in both the effective 
range expansion model ([T^ and the two-channel model. 



2.2. Model Comparison 

Before we consider the recombination rates and the binding energies in the different 
models, we first make a comparison of the two-body potential models in terms of 
their predictions for the i/^ = X + i coefhcient that provides the effective potential for 
the three-body system in the hyperradial equation ([5]). We compare the models by 
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explicitly plotting their associated eigenvalues in figure [TJ At large p all models have 
the same asymptotic value — — 2p^/a^, since we have a two-body bound state for 
a > that dominates the large distance behaviour through the configuration where 
two particle are bound and have small interparticle distance and the third particle is 
far away. This is more clearly illustrated in the inset of figure [T] where the horizontal 
axis extends up to p/atg = 1000. For intermediate distances p > 2|i?| the finite-range 
models show surprisingly similar forms given their quite different formalism. What 
is particularly important to notice is the fact that an inner pocket develops in the 
two-channel model, but at the same time a barrier with respect to the zero-range 
model is also seen. For p < \R\ the effective range expansion model goes to zero, and 
thus effectively becomes regularized. There is no need for a three-body parameter (in 
the form of a high-momentum cut-off in momentum-space or a short-distance cut-off 
in coordinate space). This is not the case for the two-channel model we use here. It 
requires a cut-off at small distance. 



3. Finite-range effects in the recombination rate 

We now proceed to consider three-body observables starting with the recombination 
rate on the positive a side of the Feshbach resonance. On this side of the resonance 
the recombination takes place by transition of two of the three particles into the 
channel with a bound two-body dimer with the universal binding energy proportional 
to — a~^. On the a < side where there is no bound dimer, the decay goes directly 
into some strongly bound two-body state of the atom-atom potential and depends on 
the short-range details. This latter case will not be considered here. 

The recombination rate for a > is calculated using the semi-classical WKB 



Finite Range Effects for Three Bosons 



8 



method of hidden crossing theory where the recombination coefficient a is 

a = 8(2^)23V3— hm ^ , 



(22) 



with the wave-number k defined hy E ~ 1i?k^ /2m*. The transition probability 
P{k) = |S'oi(fc)p, with 5oi the transition matrix element between adiabatic channels 
and 1, is [IH| 



P{k) = 4e-2Ssin2 A , 



(23) 
(24) 



where the integral is taken along a contour c in the complex p-plane connecting 
the adiabatic channel corresponding to three free particles, n = 1, to the channel 
describing a dimer and a free particle, n — ^. Additional details can be found in 
reference [3S]. Note that since i'o(O) < for the zero-range and two-channel models 
there is no classical turning point in the open (n = 1) channel, thus we employ a 
regularization cut-off in order to avoid a divergent integral. 
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Figure 2. Recombination coefficient a \22\ for the zero-range model, tiie effective 
range expansion model and the two-channel model with i?o = — Scibg (dotted) and 
_Ro = — lOajjg (full). The inset shows a closer look at the minimum near a\. Note 
that the cut-off is such that all models reproduce the minimum at a| . This allows 
us to study the effect at the other minimum of the different models. 



The recombination coefficients for different values of the effective ranges and 
different models are shown in figure[21 The cut-offs were chosen such that the minimum 
at is the same for all models. We can thus compare the models at the other 
minimum. For the zero-range model, the ratio of a\ to a\ is 22.7, showing that this 



Finite Range Effects for Three Bosons 



9 



calculation scheme agrees with the universal result. For the other models this ratio 
is reduced, the minimum at moves towards higher a. In order to make this more 
clear we plot the ratios of the minima positions as function of the effective range, 
R, in figure [3] We see that the two-channel and range expansion model give similar 
qualitative predictions but there are overall quantitative differences. We cannot extend 
the curves in figure[3]all the way to i?o = due to numerical difficulties, but the trends 
should be clear. What we also see is that the scale factor reduces quite drastically 
at large i?o for both models. This corresponds to narrow Feshbach resonances, where 
there are currently not enough experimental data to make a detailed comparison. 
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Figure 3. The ratio aj/a* for the zero-range model, the range expansion model 
and the two-channel model as function of iio/fitg- The zero-range model has a 
ratio of 22.7. This factor is reduced when the range is increased to more negative 
values. 



3.1. Cut-off effects on binding energy 

As we noted already above, we need to put a cut-off on the WKB integral in order 
to render it finite. We therefore need to consider the behaviour of the results as we 
change this cut-off. In particular, the choice of cut-off affects the atom-dimer binding 
energy. The WKB approximation yields a simple estimate of bound state energies 
(similar to the Bohr-Sommerfeld quantization rule) 



dp^2E^-!^=.[n-\^ , (25) 

where n — 1,2,... indicate the ground state, first excited, . . . etc. In the universal 
limit a — >■ cx) where vo{p) = isp, is given by 

ii;„«-^exp + -A ' (26) 

Pcut V so 2so / 
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clearly showing the geometric scaling -E„+i — e^^'"^^°En ~ i?„/22.7^ and the Thomas 
effect E„ — > — oo for pcut 0. Solving (|25p numerically for finite scattering length 
gives atom-dimer binding energies for the three models at e.g., R = —10, a = 500, n = 
2 

Zero-range. = -0.002060 , 
Two-channel. :£'b^^ = -0.002441 , 
Efl[ective range. li;^^ = -0.002062 , 

which shows that the chosen cut-offs give similar bound state energies for the different 
models. The effect of this choice of cut-off is thus under control. 




4. Bound trimers 

In order to further study the three-body physics and its dependence on two-body 
interaction models, we now consider the three-body bound state spectrum when finite 
range effects are included. When the scattering length a is large, yields the same 
binding energy values as the radial equation ([S]). However, for the following discussion 
we will solve the radial equation numerically ([5]) to obtain accurate results also when 
the binding energy is close to zero. We solve the radial equation ([5]) for f{p) such 
that /(jOmax) = for some large pmax- Also, the regularization cutoff pc is used in the 
boundary condition f{pc) = 0. 

For the three models in question we have calculated E{a) for the three lowest 
trimers, the result is shown in figure HI For positive a the dashed line indicates the 
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atom-dimer threshold which is given by the diner binding energy — 1 /a^ for the zero- 
range model and by for the effective range expansion model. For the two-channel 
model this can be calculated only numerically, yet it agrees surprisingly well with the 
analytical formula for the effective range expansion model ([15]). The cut-offs (on the 
coordinate-space hyperspherical potential) were chosen such that the second trimer 
(n = 2) energies coincide at \a\ = oo. The three spectra for n — 3 are virtually 
identical. This is reasonable since the binding energy is very small for n — 3 and the 
state is almost completely insensitive to finite-range effects. However, for the ground 
state n = 1 a clear distinction between the models appears. At \a\ = oo the finite- 
range models give practically the same trimer energy, a factor of ~ 25.3^ times higher 
than the n — 2 state. In comparison the zero-range model trimer energy is only a 
factor of 22.7^ times higher. 

We have already seen that the ratio of a-values corresponding to recombination 
minima in figure [5] differs from the universal value 22.7. This is directly related to 
the change in the atom-dimer threshold shown in figure 2] since the ratio between a 
corresponding to trimer energy termination points on the threshold line is no longer 
22.7 for the models incorporating effective range. 
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-1.9 
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Figure 5. The product K*a^~^ as a function of effective range, Hq, where 
—h?K*^/2mr = The universal value K*a'~) = —1.5076 for the zero-range 

model is not correct for the lowest trimer, where the value is —1.469, independent 
of cutoff. The two-channel model curves are for different value coordinate-space 
cut-offs on the hyperradial potential. The cut-off is 0.5, 0.6 and 0.7 in units of 
afjg for the top, middle and bottom blue curves. For the effective range expansion 
model the dependency on cutoff is insignificant. Note that the two-channel model 
only works for Rq < 0. 

For negative a the value of a^"^ indicates the threshold scattering length for 
creation of the lowest Efimov trimer, as indicated in figure |4l When written in units 
of TvdWi this quantity is the subject of much recent discussion since it seems to have a 
universal value of a^^-' ^ — 9rvdw for different cold atomic systems [15]. Here we want 
to address finite range effects on the value of a'-"-' for the two models. Some other 
recent works that address finite range effects on this threshold value can be found in 
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reference [S^ and P^ . 

Our results within the different models for a^~) as a function of Ro/abg are shown 
in figure [SI Most notable is the lowering of a^"-* for the finite-range models compared 
to the zero-range model. This is partly due to the lower binding energy at 
\a\ = oo. The product o'-^-'k* (where —h^K*^/2mr — e!^^) is universal in the zero- 
range model with the value -1.5076 [33]. Thus increasing |£^;^°°-'| will reduce la'"-*]. 
This effect is, however, not enough to account for the deviation from the zero-range 
result. The product is further reduced for decreasing Rq indicating a lower value of 
a(-). 

Both finite-range models show the same trend. However, the value of the change is 
different for the two models when —Rq gets sufficiently large. The plot also shows that 
for the two-channel model the effect depends on cut-off. The cut-off, chosen for the 
purpose of illustration but with reasonable values, is 0.5, 0.6 and 0.7 for the top, middle 
and bottom blue curves in figure [5l The effective range expansion model shows only a 
very small dependence on the cut-off (the three green curves in figure[5|show very little 
deviation). The opposite behaviour, i.e. |a^~^K*| gets larger for larger effective range, 
is found for Rq > 0. This can be connected to the use of finite range potentials with 
positive effective range. To obtain results for this case we use the boundary condition 
including the shape parameter (|16p in the effective range expansion model in order to 
extend it to positive effective range. Indeed we find that the change is now opposite, 
see figure [HI This is similar to the product of a'^^^K* found in reference [H]. We 
notice, however, that reference |22| find that the trimer binding energies for \a\ = oo 
get smaller and the value of |a'~^ | gets larger for large effective range. This is not seen 
in our calculations. 



+ 
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Figure 6. Ratio of trimer bound state energies at resonance as a function of 
effective range for several different cut-offs. The solid lines show the ratio of 
energies for the first and second state, the dashed lines for the second and third 
state. Red, yellow and blue curves (top three in legend) are for the two channel 
model while green, cyan and magenta curves (bottom three in legend) are for the 
range expansion model. 
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In order to better understand the behaviour of the trimer energies and a^~\ 
figure [S] shows the ratio of trimer bound state energies on resonance (|a| — >■ oo) for 
the two-channel and range expansion models for three different cut-offs, chosen for 
the purpose of illustration. The most noticeable feature is the rise and fall of the 
ratio of the ground state energy E^^^ to the first excited state energy for the two- 
channel model for negative i?o- This non-monotonous behaviour can be understood if 
one assumes a three-body wave function that lives at large hyperradii, p. When the 
effective range is decreased from zero {\R\ increases) a barrier with respect to the pure 
zero-range model initially decreases the binding energy. This can be clearly seen in 
figure [TJ As we increase the effective range further the wave function will leak into the 
attractive pocket at small p, which will again increase the binding compared to the 
pure zero-range result. This effect is strong for the ratio of the two lowest trimers but 
becomes weak for the ratio of the two highest trimers. This is understandable since 
the least bound trimer resides at very large hyperradii and is largely insensitive to the 
changes in the hyperradial potential at small p due to the effective range correction. 

With this interpretation the behaviour can now be better understood. In figure [7] 
we plot the solution to (jlj which determines the effective three-body hyperradial 
potential for the different model and for different signs of the effective range. The 
finite range corrections for the two-channel model and the effective range models are 
very different as one has the inner pocket and the other does not. This is the origin 
of the differences seen in figure ID Initially we see a repulsive effect compared to the 
pure zero-range model that drives the ratio upwards as all states are pushed out to 
large distance. However, for the two-channel model we eventually feel the presence 
of the inner pocket and states will leak into smaller hyperradii where the ratio of 
energies is in turn driven down. For the positive effective range case we see from 
figure[7]that the range expansion model will now have a pocket at small distance. The 
trimer states can leak into this pocket and the ratio of energies will again go down. 
The important point is that in the models we have presented here there is no reason 
to expect monotonic behaviour since the effective three-body hyperradial potentials 
have non-trivial structure with repulsive and attractive parts in comparison to the 
pure zero-range model with no effective range corrections. 

5. Conclusion 

We investigate finite-range effects in three-body recombination rates in cold atomic 
gases near Feshbach resonances as well as finite-range effects in the trimer bound 
state energy spectrum. We use two models which include the finite-range effects 
and compare their results with a contact interaction (i.e. zero-range) model. The 
first model is the range expansion model which is a straightforward extension of the 
zero-range contact interaction model. Here the effective range is included directly in 
the boundary condition on the three-body wave function following the effective range 
expansion of standard scattering theory. Variation of the scattering length through 
the Feshbach resonance is done phenomenologically as in the zero-range model. This 
model can also be used for positive effective range calculations. The second model is a 
two-channel contact interaction model which naturally includes both the finite effective 
range and the variation of the scattering length through the Feshbach resonance. 

We show that with these well-tested two-body interaction models the three-body 
physics can display complicated non-monotonic behaviour as the effective range is 
varied. In particular, we find that the geometric scaling factor of 22.7 for equal mass 
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Figure 7. The eigenvalue solutions to (|4]l as a function of hyperradius, p for tlie 
zero-range , two-channel model and range expansion model for a representative 
case with |_Ro| = 5aj,g. The lowest adiabatic potentials multiplied by as 
functions of hyperradius, p, for the zero range one-channel model and three 
different effective ranges for the two-channel model. 



particles changes when including range corrections, and that it can become both larger 
and smaller than 22.7 depending on the magnitude and sign of the effective range. 

In the current setup this can be understood based on the functional form of the 
effective hyperradial potential. On resonance where the scattering length diverges, 
the lowest trimer bound state has the strongest dependency on effective range since it 
lives at small hyperradius, whereas the excited states live at much larger hyperradii 
and the effective range contribution is much less profound. The adiabatic potential 
of the range expansion model is raised and lowered relative to the zero-range model 
potentials when the effective range is negative and positive respectively. This leads to 
bound states being less or more bound respectively. For the two-channel model the 
effective range is always negative which can only be achieved by using potentials with 
an outer barrier. The hyperradial potential reflects this fact and develops a pocket at 
small hyperradii that the lowest states will eventually leak into. This feature is similar 
to the range expansion model but for the case of positive effective range. 

Our results demonstrate that effective range corrections within the framework 
of contact model potentials can lead to non-trivial behaviour of trimer energies, 
thresholds and interference features in recombination rates. Effective range corrections 
are expected to be important for the case of narrow Feshbach resonances [10]. The 
experimental data on Efimov states for narrow resonance systems is sparse and more 
measurements are needed in order to fully discriminate between different models that 
include finite range corrections. However, what we can conclude is that care must be 
taken when a particular two-body scattering model is used for the trimer states that 
have the largest binding energies in a universal setup, i.e. for the lowest states that 
have binding energies related to the background short-range scales. For higher lying 
trimcrs it is less important since the states are largely insensitive to the short-distance 
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behavior of the effective three-body potential. 
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